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ABSTRACT 

In this work we study how the cosmological parameter, the Hubble constant H , 
can be constrained by observation of very high energy (VHE) 7-rays at the TeV scale. 
The VHE 7-rays experience attenuation by background radiation field through e + e~ 
pair production during the propagation in the intergalactic space. This effect is pro- 
portional to the distance that the VHE 7-rays go through. Therefore the absorption of 
TeV 7-rays can be taken as cosmological distance indicator to constrain the cosmolog- 
ical parameters. Two blazars Mrk 501 and 1ES 1101-232, which have relatively good 
spectra measurements by the atmospheric Cerenkov telescope, are studied to measure 
H . The mechanism measuring the Hubble constant adopted here is very different 
from the previous methods such as the observations of type la supernovae and the cos- 
mic microwave background. However, at 2a level, our result is consistent with which 
given by other methods. 

Subject headings: galaxies: distances and redshifts — BL Lacertae objects: individual 
(Mrk 501, 1ES 1101-232) — cosmological parameters — gamma-rays: general 



Introduction 



The modern cosmology achieves great progress in recent years. A concordance ACDM 
cosmology has been built thanks to th e precise observations of the " distance indicators" type la 



supernovae (SNe la, IRiess et all Il998l 12004c IPer lmutte r et al 



mic microwave background (CMB, Ide Bernardis et al.ll200ol: Ispereel et al. 
other cosmologica l probes such as the large scale struct ures (ITegmark et al. 



1999 1) and the an isotropy of cos 



2003). There are also 



2004|) . galaxy clusters 



(| Allen et al.l 120081). observationa l Hubble parameters (|Simon et al.l |2005|) and the weak gravita- 
tional lensing (|Munshi et al.ll2008|) further supporting this scenario. Different methods are roughly 
consistent with each other within the observation uncertainties. It is very important to develop 
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additional complementary observational evidence to test this model and measure the cosmological 
parameters. 

It has been pointed out that the observations of VHE 7-rays at the energy scale of hundred 
GeV to TeV scale are possible to provide another independent constraint on the cosmological 
parameters (|Salamon et al.l 1 19940 . Thanks to the rapid technical development of VHE 7-ray de- 
tection, especially the atmospheric Cerenkov telescopes, great progress of VHE 7-ray astronomy 
is achieved in recent years and a large number of VHE 7-ray sources are detected. Even 7-ray 
sources at cosmological distances, such as Mrk 501 at z = 0.034 and 1ES 1101-232 at z = 0.186 
are observed. More importantly the spectra of these sources have been measured with relatively 
high precision, which provide us the possibility to untangle the effect of attenuation when 7-rays 
propagate in the intergalactic space. 

The attenuation of VHE 7-rays is induced by the electron -positron pair production 7 + j hk — » 
g + + g~ durin g its propagation in the background radiation field (|Nikishovll 19621 : iGould & Schredei 
19671 1 1966|) . This process is actually complex. The observed spectra of extragalactic sources are 
related with several issues: the intrinsic spectra at sources, the cross section of 77 interaction, the 
intensity of the cosmic infrared background (CIB), and the physical distance the VHE 7-rays cross. 
Even before the first detection of the VHE 7-rays from distant e xtragalactic sources , the perspective 
to explore the CIB using the attenuation effect was proposed (|Stecker et al.lll992|) . The discovery 
of the first extragalactic VHE 7-ray source, an active galactic nuclei (AGN) Mrk 421, was per- 
formed by Whipple in 1992 (|Punch et al.lll992h . Till now more than 20 extragalactic sources, 
most of which are AGNs, are discovered by ground-based observatorie^. The observations of 
these sources provide us valuable information in understanding the 7- ray production mechanism 
and give useful implication or constraint on the CIB intensity (e. g., Coppi & Aharoniar3ll999l : 
Krawczynski et alj|2000c iRenault et al.ll200ll : lAharonian et al.ll2006|) . On the other hand, once the 
primary spectra and the CIB intensity are specified, the distance-redshift relation (accordingly the 
cosmological model parameters) of the sources can be derived from the absorption effect. 

In this work we try to constrain the Hubble constant from the absorption effect of distant VHE 
7 sources by the CIB. By a global fitting to the observational spectra of two TeV blazars, Mrk 501 
and 1ES 1101-232, we get the Hubble constant with larger errors compared with other methods. 
In our w ork the ACDM unive rse with matter component Vl M = 0.28 and dark energy Qa = 0.72 is 
adopted (IKomatsu et al.ll2008|) . We find that the best-fitting to the data of the two sources intend 
to give similar Hubble constant, although they have very different intrinsic spectra and redshifts. 
This is very encouraging that the attenuatio n may indeed give implications on the cosmological 
parameters. We noticed in a previous work lBarrau et al.1 (|2008f) adopt the similar effect to derive 
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the lower limit of the Hubble constant from observation of Mrk 501. In their work, the direct 
measurements of CIB intensity was adopted and the intrinsic spectrum of the source was required 
to be concave. 

The outline of this paper is as follows. Sec. 2 describes the absorption of VHE 7 photons 
by CIB. In Sec. 3 we present an introduction to the observations of the two TeV blazars. The 
implication on Hubble constant is given in Sec. 4. Finally we give conclusion and some discussion 
in Sec. 5. 



2. Absorption of TeV 7-rays in CIB 



The fundamental process of the VHE 7-ray absorption is due to electron/positron pair pro- 
duction 7+ 7m — > e + + e~. The threshold energy of the pair production is m 2 e /e, with e the energy of 
the background radiation. For the CMB photon e ~ 10" 3 eV, this absorption takes place for 7-rays 
with energy E > 1 PeV. While for the TeV scale 7-rays that the current experiments can probe, the 
responsible soft photon is in the infrared band, i.e., CIB with e ~ 1 eV (A ~ 1 /mi). An approximate 
relation between energies of attenuated VHE 7-rays and the CIB photons is 



A 



1/im 



1.2- 



lTeV 



(1) 



The observed VHE 7-ray spectrum after attenuation is given by 

^obs ^ -^int 5 



(2) 



where r is the optical depth and F im is the intrinsi c spectrum at the sourc e. For the CIB with 
number density n(e), the optical depth r is given as (IGould & Schredenll967|) 



T{E) 



1 — cos 

dl I dcos6> / den(e)cr(£',e,cos6 l ), 



(3) 



where dl = cdt = 



,, , dz = is the differential path traversed by the VHE 7-rays, 9 is the 

(l+z)V - 28 a+z) 3 +0-72 

angle between the momenta of VHE 7-ray and CIB photon. The cross section of pair production 
is 



&{E , 6, cos 0) = a t 

2s 
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Fig. 1. — Attenuation factor of VHE 7-rays for sources Mrk 501 and 1ES 1101-232. The three 
curves for each source correspond to the "nominal" (middle curve), 25% higher (lower curve) and 
25% lower (upper curve) CIB respectively. The Hubble constant in the calculation is adopted as 
h = 0J. 



with oj = 6.65 x 10~ 25 cm 2 the Thomson cross section, s = 2Ee(l -cos#)(l +z) 2 the center of 
momentum system (CMS) energy square, E e = and p e = y/E* — m* the CMS energy and 

momentum of electrons. 

From Eq. © we can see that the intensity of CIB is crucial in determining the effect of 
attenuation. The CIB is generated by stars and absorption/re-emission of star light by dust in 
galaxies. The status of measurements and models of CIB can be found in the review paper by 
Hauser & Dwekl (|200lh . Because of the contamination of foreground from the solar system and 
the Galaxy, the deter mination o f CIB has relative large uncertainty. Here we adopt the "nominal" 
model prediction of lAharonianl (|200ll) (curve 1 of Fig. 1) whic h can give a good description of 
the measurements. Two other models (curves 2 and 3 in Fig. 1 of lAharonianl (120011) ) are regarded 
as the lower and upper limits of CIB intensity. The differences between these model predictions 
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can vary from several tens percent to several times at different energies. To simplify the uncer- 
tai nties of CIB in ou r analysis we take the uncertainty of ±25% relative to the "nominal" model 
of lAharonianl (1200 1|) to repres ent the upper and lower limits. The CIB is denoted as n(e) =An(e), 
where n(e) represents the best ( Aharonianll200 1 ) model of CIB, A = 1 ± 25% is a normalization fac- 
tor to represent the uncertainties, which is energy independent. This form of uncertainties greatly 
simplifies the process of global fitting. 

The comoving density of CIB is adopted to be constant witho ut redshift evolution, w hich is 
shown to be of little influence for the sources with redshift z < 0.2 (|Aharonian et al]|2006l) . Using 
this CIB field, we calculate the attenuation factor e~ T of VHE 7-rays according to Eq.® for sources 
Mrk 501 (z = 0.034) and 1ES 1101-232 (z = 0.186), as shown in Fig. [B It can be seen from this 
figure that the absorption for 7-rays increases rapidly for energies > 10 TeV. It also shows that 
absorption of the nearby source Mrk 501 at energies ~ 20 TeV is comp arable with the effect to 
the distant source 1ES 1101-232 at energies ~TeV (IBarrau et al.1 12008). Therefore 7-rays with 



high energy (E > , 10 TeV) or high red shift (z > 0. 1) will be very effective to study the attenuation 
process and CIB (lPrimack et al.lll999l) . 



3. TeV 7-ray Observations of the sources: Mrk 501 and 1ES 1101-232 



Mrk 501 is a nearby (with redshift z = 0.034) BL Lac type blazar, which is a kind of radio-loud 
AGN with relativistic jet being aligned along the line of sig ht. The first detect ion of VHE 7-ray 



emission from Mrk 501 was performed by Whipple in 1995 (Ou inn et al. 



flares , Mrk 501 was observed by several experiment groups dC atanese et al 



19981 : lAharonian et all 1 19971 : bjannati-Atai et all 1 19991 : lHavashida et al. 



996). During the 1997 



1997; Samuelson et al. 



1998|). The energy up to 



~ 20 TeV observed from Mrk 501 makes it a good c andidate to study the absorption effect of 
high energy 7-rays in CIB (ICoppi & Aharonianlll999|) and the possible Lorentz violation effect 
(|Protheroe & Meyerll2000|) . I n the current study, we use the reanalyzed HEGRA data in 1997 with 
improved energy resolution (lAharonian et al.l 1200 1|) . The observational spectrum is shown in the 
left panel of Fig. [2] 

The other source investigated here, 1ES 1 101-232, is a distant blaza r with redshift z = 0. 1 86. 
TeV observation of 1ES 1 101-232 was performed by H.E.S.S. in 2005 (lAharonian et allbood) . It 
is shown that even though the measured maximum energy only reaches ~ 3 TeV, 1ES 1 101-232 is 
still very effective to probe th e intensity of CIB in the propagation path VHE 7-rays go through due 
to its large distance from us (lAharonian et al.ll2006|) . The observational spectrum of 1ES 1 101-232 
is shown in the right panel of Fig. [2] 



Actually these observations show harder spectra than expected. Considering the effects of 
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Fig. 2.— Left: Energy spectrum of Mrk 501 by HEGRA in 1997 (lAharonian et all 1200 lh . Also 
shown are the absorption-corrected spectra using the usually adopted Hubble constant h = 0.7 
(upper most points) and the best-fitting Hubble constant h = 0.98 (medium points, see below §4). 
The two lines are the best-fitting power law intrinsic spectrum (thick-black, 1^ = 2.62) and the 
absorbed one (thin-red ). Right: same as the left panel but for 1ES 1101-232 by H.E.S.S. in 2005 
(lAharonian et al.ll2006|) . The best-fitting parameters are h = 1.09 and r int = 0.93 instead. 



absorption of the VHE 7-rays from the normal CIB density and distance-redshift relation, the ob- 
served spectra means unnaturally hard intrinsic spectra at the sources. It can be clearly seen from 
the upper most points in Fig. [2l which are calculated using the usually adopted CIB and h = 0.7. 
For 1ES 1101-232, the corrected spectrum is ~ E~° A , which seems extremely har d when compar- 
ing with the expected one from shock acceleration (IMalkov & O'C Druryll200l|) . The corrected 
spectrum of Mrk 501 can not be fitted with a single power law, however, it is also shown in Fig. |2] 
that the spectrum of high energy part is very flat. 

This possible anomaly has led to quite a few discussions about possible new physics. It was 
proposed that the axion-7 oscillation when the VHE 7-rays propa gate in the interg alactic magnetic 
field makes the universe more transparent than naively expected (|simet et"aL 2008^) . When 7-rays 
oscillate into axions they will not be absorbed by the CIB and keep the primary spectra unchanged 
(|Hooper & Serpicoll2007l) . ] t was also suggested that the p ossible Lorentz violation may be respon- 
sible for the hard 7-ray spectra ( Protheroe & Meyer! |2000|) . In this scenario the threshold energy 
of the interaction moves to higher energy and the absorption effect at the observed energy scales 
becomes weaker. Possible ex plan ations of the hard spectr a within astrophysics are also discussed 
(lAharonian et alJbood 12008b . In lAharonian etall (120061) the authors pointed out that if the CIB 
intensity is about half of the locall y measured values the observed 7-ray spectra can be naturally 
explained. lAharonian et al.l d2008|) also suggested some special mechanism to produce very hard 
intrinsic spectra at the sources. 
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Since all these explanations are based on a standard cosmological model, we are considering 
that if these observations have implications on the cosmological model itself. In the following we 
will show that the best-fitting to the 7-ray data favors a larger value of the Hubble constant. In 
spite of the large uncertainties, our result is consistent with previous cosmological measurements 
at 2a level. 



4. Implication on the Hubble constant 

The 7-ray spectra from astrophysical sources are usually very well described by power law 
functions, which originated from the shock wave acceleration at the sources. Assuming the power 
law spectral index ri nt at source we get the observed spectrum F obs oc E~ Tmt e" Ar °/ /l according to Eq. 
©, where r is the optical depth with A = 1 and h = 1. Using the observational data we can fit 
the parameters r int , A and h. It should be noted that parameters A and h are strongly coupled with 
each other, so it is unable to determined them simultaneously from the attenuation of VHE 7-ray 
spectra. We firstly fix A = 1, and fit the parameters r int and Hubble constant h. Then we will take 
the uncertainty of A into account. 





2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 



0.5 1 



Fig. 3. — 1, 2 and 3a confidence regions (from inner to outer) of parameters r; nt and h for A = 1. 
Left: for Mrk 501; right: for 1ES 1 101-232. The cross in each panel is the best-fitting value. 



The confidence regions in the h-T plane are shown in Fig. |3] The best-fitting values and la 
errors of the parameters are compiled in Tabled! It is shown that for both sources, the best-fitting 
Hubble constant h is close to 1, which i s larger than the results from other cos mological mea 
surements h « 0.7, such as from SNe la (lAstier et al.ll2006|) and CMB anisotropy (IKomatsu et al. 



2008|). A larger h implies that smaller abso rption is favored by the observations. Similar re sults 
are also found by the previous studies (e.g.. |Protheroe & MeyerM2000t lAharonian et alJ l2006). We 
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also notice in Fig. [3] that since the statistical errors of Mrk 501 are much smaller than that of 1ES 
1 101-232 it also gives much better constraints on the parameters. 

The best-fitting intrinsic spectrum is T int = 2.62 for Mrk 501 and 0.93 for 1ES 1101-232. 
It seems that the intrinsic spectrum for 1ES 1101-232 is still too hard. Generally the intrinsic 
spectrum of blaz ars for both hadronic and leptonic scenarios from shock acceleration is expected 
to be r int > 1.5 JMalkov & O'C Druryll200ll : lAharonian et alJl200fih . If we apply a limit r int > 1.5, 
we find that the Hubble constant h > 1.2 at 68% confidence level from 1ES 1101-232. It should 



be no ted that scenario with very hard 7-ray spectrum is also proposed recently (|Aharonian et al 
20081) . 



It is also shown in Fig. [3]that there is degeneracy between the parameters T^t and h, especially 
for the source 1ES 1 101-232. This is because a harder r int means a stronger absorption, and leads 
to a smaller h (or a larger A). 

The reconstructed source spectra using the best-fitting parameters are shown by the medium 
points in Fig. [2] The thick lines in this figure represent the best-fitting power law intrinsic spectra. 
We can see that the reconstructed source spectra are well consistent with power law functions. 




2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 




Fig. 4. — Same as Fig. [3]but with a prior A = 1 ± 0.25. 

Furthermore, we take the uncertainty of the CIB into account by employing a prior A = 1 ± 
0.25 when doing the global fitting. The fitting results of r int and h are shown in Fig. HI For 
Mrk 501, the best-fitting values are almost the same as the case A = 1 (Fig. [3]), but the contours 
become larger after including the uncertainty of CIB. The fitted Hubble constant with la range is 
h = 1. 01^ 4Q. While for 1ES 1 101-232, the best-fitting values of the parameters differ significantly 
from the case A = 1 and have larger uncertainties, as shown in the right panel of Fig. 0] 

Finally, we combine data of the two sources to fit the Hubble constant h. We show the x 2 
values as functions of h taking the prior A = 1 and A = 1 ± 0.25 respectively in Fig. [5J The lines 
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Fig. 5. — x 2 as functions of h for the combined fitting using both Mrk 501 and 1ES 1 101-232. Left: 
A = 1; right: A = 1 ±0.25. The two horizon lines in each panel are x 2 = Xmm + 1 an ^ X 2 = Xmin + 4> 
which represent the la and 2a uncertainties of h. The shaded regions show the 2a results on the 
Hubble constant h = .72 ± 0. 16 (wider one) from the Hubble Space Telescope (HST) Key Project 
( Freedman et alJl200l|) and h = 0.701 ± 0.026 (narrower one) from the recent observations of CMB, 
SNe la and baryon acoustic oscillation (|Komatsu et al.l l2008). 



X 2 = Xmm + 1 (X 2 = Xmin + 4) is plotted to show the la (2a) range of parameter h. We find that 
h = LOO!^ for A = l,h = 1.05^9 for A = 1 ±0.25 at la level respectively. We can see that after 
combining data of the two sources the best value of h is not sensitive to the uncertainties of CIB, 
although the error bar of h becomes larger. 



5. Conclusion and discussion 

In this work we constrain the cosmological parameters, especially the Hubble constant H , by 
observations of extragalactic VHE 7-ray sources at cosmological distances. The VHE 7-rays ex- 
perience attenuation by background radiation field through e + e~ pair production. This attenuation 
is proportional to the distance that the VHE 7-rays go through. Therefore the absorption of VHE 
7-rays can be used to determine the distance of VHE 7-ray sources, accordingly to get constraints 
on the cosmological parameters. 

By fitting the spectra of two blazars Mrk 501 and 1ES 1 101-232 we get the best-fitting Hubble 
constant is H ~ 100 km s" 1 Mpc" 1 . A large Hub ble constant implies that the absorption of VHE 7 - 
rays is not as significant as we usually expected JProtheroe & Meyeill2000 ; Aharonian et al. 2_006). 



Since the observations of VHE 7-rays and CIB are still rough, the errors of the fitting parameters 
are also very large. The mechanism constraining the Hubble constant adopted here is very different 
from previous methods, however, our results are consistent with the recent combined analysis of 
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CMB , SNe la and baryon acoustic oscillation data on Hubble constant h = 0.701 ±0.026 at 2a 
level (komatsu et alJbood) . 



In fact, for each single method measuring the Hubble constant there are relatively large un- 
certainties, including both the statistical and the systematic ones. The HST Key Project mea- 
sured the Hubble con stant from several secondary distance indicators using Cepheid as calibration 
JFreedman et alJboOlk . They gave the results that h = 0.71 ± 0.06 for SNe la, h = 0.71 ± 0.08 for 
the Tully-Fisher relation of spiral galaxies, h = 0.70 ±0.08 for the surface brightness fluctuations 
of galaxies, h = 0.72 ± 0. 1 1 for Type II supernovae and h = 0.82 ± 0. 1 1 for the fundamental plane 
method of elliptical galaxies. The combined result of HST Key Project suggested h = 0.72 ±0.08. 



Whil e after the CMB observations are involved, the result is greatly improved (IKomatsu et al 



2008J), as shown in Fig. [5J It shows that the cross check and combination of different methods are 
very helpful to find the right answer and improve the accuracy. 

Before concluding we would like to briefly comment the simple assumptions adopted in our 
work. We assume the intrinsic spectrum to be a single power law. In fact the spectral energy dis- 
tribution of many AGNs can be well described by the so-called synchrotron self-Compton (SSC) 
model. In the SSC scenario, the VHE 7-rays are produced through the inverse Compton (IC) scat- 
terings between the electrons and synchrotron photons generated themselves. The VHE spectrum 
generall y has an "IC peak" originat es from the transition from Thomson regime to Klein-Nishima 
regime (|Blumenthal & Gouldlll970|) . However, if the energy range is narrow, e.g., within a decade 
of energy, it can be described approximately by a power law. 

The CIB model and its uncertainties are also too simplified. Because of the contamination 
of foreground radiation, it is usually difficult to get reliable CIB from measurements. The galaxy 



evolu tion models to predict C IB also have large uncertainties (IHauser & Dwekll200ll : IStecker et al. 



20061) . In lStecker et al.l(|2006l) the results between t he "baseline" and "fas t" evolution models differ 
by ab out 20% ~ 40%. While the results given by IPrimack et al] (|2005|) differ from IStecker et al. 
(|2006|) by a factor of 2 at some wavelengths. 



Anyway we think the present work is only a prototype of such studies, since the observation 
of VHE 7 rays from extragalactic sources is achieved only in the recent years. This field is actually 
immature and at its early stage. With the next generation of space and ground-based instruments 
more extragalactic sources will be observed with high precision. Especially the space observatory 
Ferm^l (energy range from MeV to hundred GeV) can explore 7-ray sources to redsh ift z ~ 1, due 
to an estimate of the "7-ray horizon" log(z) ~ l-O^logCE/lGeV) (lHartmannll2007l) . The "7-ray 
horizon" is the redshift corresponding to absorption depth r«l for energy E. With larger sample 
of data, higher precision of spectra and higher redshift sources from Fermi we can even explore 



2 See the homepage of Fermi, http://www-glast.stanford.edu 
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more cosmological parameters besides the Hubble constant, such as the cosmological component 
or the equation of state of dark energy. The development of ground-based instruments also aims 
to lower threshold energy and improve sensitivities. We anticipate the field of VHE 7-ray will 
develop quickly in the near future. Our work shows that the VHE 7-rays may become an more 
important field, not only to astrophysics but also to the cosmology. 

This work is supported by the NSF of China under the grant Nos. 10575111, 10773011 and 
supported in part by the Chinese Academy of Sciences under the grant No. KJCX3-SYW-N2. 
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Table 1 . Fitting results for parameters 



Source A = 1 A = l±0.25 



Mrk 501 

1ES 1101-232' 

Combined 



2 62+ - 11 
0.93 



1.09 



2 fS2 +(U2 
2.15 



ft 

1 01 + - 53 

2.80 
1 05+ - 35 



a The fitting errors of parameters for this source are very 
large that not shown here. 



